A model of hydrodynamic interaction between swimming bacteria 
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We study the dynamics and interaction of two swimming bacteria, modeled by self-propelled 
dumbbell-type structures. We focus on alignment dynamics of a coplanar pair of elongated swim- 
mers, which propel themselves either by "pushing" or "pulling" both in three- and quasi-two- 
dimensional geometries of space. We derive asymptotic expressions for the dynamics of the pair, 
which, complemented by numerical experiments, indicate that the tendency of bacteria to swim in 
or swim off depends strongly on the position of the propulsion force. In particular, we observe that 
positioning of the effective propulsion force inside the dumbbell results in qualitative agreement 
with the dynamics observed in experiments, such as mutual alignment of converging bacteria. 

PACS numbers: 87.16.-b, 05.65. -|-b, 87.17.Jj 

I. INTRODUCTION 

Modeling of bacterial suspensions and, more generally, of suspensions of active microparticles has recently become 
an increasingly active area of research. One of the motivating factors behind this trend is the study of the dynamics 
of large populations of aquatic single-cellular [TJ [SJ [31 SJ |S] and multicellular organisms [51|7]. In particular, there is 
significant interest in understanding the mechanism of formation of coherent structures on a scale much larger than 
individual microorganisms in suspension (see, e.g., [5l[8l[9j). Recently, bacterial suspensions have also emerged as a 
prototypical system for the study and engineering of novel biomaterials with unusual rheological properties [lOl [11] . 
Here the idea is to exploit the active nature of the particles in the suspension in order to generate specific effects, such 
as enhancement of transport and diffusion of tracers relative to that of the solvent [5,6, 10 . A good review of the 
motivations, experimental studies, and modeling approaches to suspensions of swimming microorganisms is contained 
in the introduction to [T^ . 

The principal organizing role in the formation of large-scale patterns (e.g., [5J [31 [5J [T3]) is believed to be played by 
hydrodynamic interactions between individual swimmers and the environment. This includes the boundary effects as 
well as the hydrodynamic interaction with other swimmers [H [8] . These effects and interactions are also believed to set 
the spatial and temporal scales of the patterns. At the same time, fundamental questions about the hydrodynamics of 
a single swimmer have been studied by many researchers over several decades (e.g., [T51ITB1 [TTI [T5] and references 
therein). Here one of the central features is the very low Reynolds number TZe of a typical microscopic swimmer 
- TZe ~ 10^"* - 10^^ [m [16] - making the governing dynamics Stokesian. Since the Stokesian dynamics is time- 
reversible, the very possibility of propulsion at low Reynolds numbers had to be clarified in general (see, e.g., [17]), 
with some of the early important contributions made by Purcell [19^ . Specific studies of the propulsion mechanism of 
fiagellates includes the work by Phan-Thien et al. [201 [H]. In particular, Ref. [2U| in detail addresses the hydrodynamic 
interactions of two nearby microswimmers. We also address the question of pairwise hydrodynamic interactions of 
swimmers. Unlike in |20j . however, our model abstracts from the method of propulsion (e.g., rotating helix, water 
jet) and is applicable to a wider class of swimmers. Also, being simpler structurally, our model allows us to perform 
simulations for a larger collection of swimmers and relatively long time. 

Studies of the fundamental interactions of small numbers (e.g., pairs) of swimming particles are also important for 
validating mean-field theories of large-scale pattern formation in active suspensions f22J . The continuum phenomeno- 
logical models proposed in such studies typically rely on a two-phase formulation of the problem: the particle phase 
interacts with the fluid phase via a postulated coupling mechanism. It should be possible, at least in principle, to 
derive or verify the proposed coupling mechanisms against the fundamental particle-particle dynamics. In this paper 
we focus on the mechanisms of alignment of a pair of elongated swimmers, with the aim of shedding light on a possible 
mechanism of large-scale ordering in dilute bacterial suspensions. 

The main emerging approaches to the modeling of microswimmers [121 [18l [231 El] typically abstract away the 
details of the actual propulsion mechanism and use simple, tractable, rigid geometries to model a swimmer. In 
general "higher-order" effects such as signaling between bacteria and chemotaxis are ignored, the emphasis being on 
the basic hydrodynamic interactions. Both [25 and [24 model swimmers as elongated bodies; the former employs 
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slender body theory for cylindrical rods, while the latter models the elongated body as a dumbbell consisting of a 
pair of balls. As our model is a modification of [21], the following sections contain a more detailed description of 
the dumbbell model. In the aforementioned studies the self-propulsion mechanism is modeled by a prescribed force, 
concentrated at a point inside the dumbbell ball, as in [2?, or distributed over a part of the surface of the body in the 
form of a specified tangential traction, as in 23,- Both of these studies ultimately rely on numerical simulations 
with the goal of studying the emergence of large-scale coherent patterns predicted by continuum theories such as 
[51 122| or observed in experiments [3^, "F. By contrast, the main tool of our work is an asymptotic analysis followed by 
straightforward numerical simulations. 

A different structural and dynamic approach is taken by Pedley and coworkers (see [12 US])- Here a basic swimmer 
is modeled as a squirming sphere, with a prescribed tangential velocity as the model of the propulsion by motile 
cilia (short hair on the surface of the cell beating in the same direction). As with the dumbbell model, a spherical 
squirmer allows using some fundamental solutions and relations (e.g., the Stokeslet solution and the Faxen relations) 
to approximate the dynamics of the swimmers. The work [12 is closer to ours in its goal of quantifying the interaction 
of a pair of swimmers, rather than a large collection of swimmers, while the model in |24j is closer to ours in the 
structural model of a swimmer. 

In our work, self-propulsion is modeled by prescribed propulsion forces (as in .2^1 vs. prescribed velocities on the 
boundary, as in [12] ; see also recent work [26] on the rheology of bacterial suspensions) , and an elongated body of 
bacterium is modeled by a dumbbell as in |24j . Our model is consistently derived from the equations of Stokesian 
fluid dynamics. In particular, we model self-propulsion by a point force, whose location can vary. This allows us to 
investigate the dependence of mutual dynamics (swim in/off) of neighboring bacteria on the position of this force 
(which roughly can be interpreted as the effect of the shape of the microorganism and the way of propulsion or 
distribution of cilia). We study the hydrodynamic interaction for well-separated swimmers. For this reason we say 
that two swimmers swim in, starting from a given mutual orientation, if at some point the distance between them 
decreases to the order of their size, (i.e., they become not well separated). Two swimmers, starting from a given mutual 
orientation, swim off if the distance between them increases to infinity without swim in happening first. We observe 
that positioning the propulsion force between the dumbbell balls results in attractive behavior of swimmers. On the 
other hand, positioning the propulsion force outside the dumbbell results in repulsive behavior. For comparison, in 
the earlier work [23j the position of the propulsion force was fixed (center of a ball in a dumbbell) . 

In this work we study the alignment of a coplanar pair of three-dimensional elongated swimmers, which propel 
themselves by "pushing" or, "pulling", mimicking a variety of self-propelled microorganisms, from sperm cells and 
bacteria to algae. We derive asymptotic expressions for the dynamics of the pair, which, complemented by numerical 
experiments, indicate that the tendency of bacteria to swim in or swim off strongly depends on the position of the 
propulsion force. In particular, we observe that positioning of the propulsion force inside the dumbbell results in the 
qualitative agreement with the dynamics observed in experiments [9] . We also observe that the dynamics of bacteria 
in a thin film (with no-slip boundary conditions on the top and bottom) is qualitatively similar to that for the whole 
space. 

One of our objectives is to develop a well-posed PDE model of an active suspension derived from first principles 
(unlike many engineering models that use ad hoc assumptions). Our proposed model is simple enough to allow for 
theoretical analysis (asymptotics) yet captures basic features observed in experimental studies. 

The paper is organized as follows. In section [llj we derive a full PDE model for the dynamics of swimmers based on 
Stokesian hydrodynamics. The well-posedness of the problem was demonstrated and can be found in Appendix [B] 

In section [1X1] we introduce an asymptotic reduction of the PDE model in the dilute limit of swimmer concentration. 
Here we show how to solve the reduced model numerically and give an analytic (asymptotic) solution for a pair of 
swimmers. Then we analyze two basic (physically interesting) configurations of swimmers based on the asymptotic 
formulas and numerical calculations. In section |IV| we make some concluding remarks and indicate areas for future 
study. Appendix [X] contains the fundamental solutions to the Stokes problem, which are used extensively in section 
|III| Appendix |C] contains some technical asymptotic formulas and calculations. Appendix [D] contains the stability 
analysis of certain configuration of bacteria. 

II. THE MODEL 

While the modeling of suspensions of passive particles is a well-established area, the mathematical study of sus- 
pensions of active swimmers is still an open area without universally accepted models that can serve as benchmarks 
for analytic and numerical studies. In this paper we propose a model that can serve as a tractable reference case 
for the mathematical questions such as existence and uniqueness of solutions. We analyze this model by considering 
its asymptotic reduction in the far-field regime and comparing its predictions with experimental results for bacterial 
suspensions. The main goal is to establish a model amenable to an analytic treatment that still captures the main 
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physical effects, such as alignment and the emergence of large-scale coherent structures; here we focus on the question 
of the pairwise dynamics of swimmers. 

We set up a PDE model for a collection of swimmers and then consider an asymptotic ODE reduction for a pair of 
swimmers that is suitable for numerical analysis. We address the question of the solvability of the PDE system and 
carry out a numerical study of the ODE model. 



A. Structure and dynamics of a single swimmer 

Structurally, we model as a dumbbell (see Fig. [T]) the elongated body of a bacterium (e.g.. Bacillus subtilis), referred 
to as the swimmer below. Similar approximation was used in [24 . It consists of two balls of equal mass m and radius 
R rigidly connected to one another a distance 2L apart. Assuming a high aspect ratio of a swimmer, we have 2L 3> R. 
The balls are denoted (head) and B^ (tail), with their centers located at x^ and x^, respectively. The unit vector 
directed from x^, to x^ is denoted by r = (x^ — Xy)/|x^ indicating the direction of the swimmer's motion, 

and the line connecting the centers of the balls is referred to as the dumbbell axis. The action of the flagellum - the 



Flagellum modeled 
by a localized force 




of the bacterium 
connected by a rigid rod 

FIG. 1: Model of a single bacterium: two balls (labeled head and tail) and the propulsion force (red ball with arrow) connected 
by a rigid rod (that does not interact with the fluid). 



bacterial propulsion apparatus ~ is represented either by a smooth volume force density F supported in a ball Bp oi 
small radius g (p <C 1), with the center Xj, located on the dumbbell axis, or by a delta function concentrated at x^. 
This force and its support will sometimes be referred to as the propulsion force. The location Xp of the propulsion 
force relative to the positions x^ and x^ of the balls in the dumbbell is defined by the parameter 

^ ^ i^p-^c)-r ^ ^^^^^ ^ ^H+x^ ^ 

For instance C — corresponds to the the center of the tail or head ball, respectively, and C = corresponds to the 
the center of the dumbbell. 

Depending on the value of swimmers are classified into pushers/pullers and inner/outer swimmers. A swimmer 
is called a pusher if ^ < and puller if C > 0. A swimmer is called outer if |(^| > 1 and inner if |C| < 1. 
The total force exerted in the ball B^ has magnitude fp = const and is directed along the axis t, that is. 



F(x)dx = -Fp = -/pT. (2) 



Bp 



Here — Fp is the force of the flagellum on the fluid, and Fp is the force of the fluid onto the flagellum, that is the force 
that propels the dumbbell. For simplicity of presentation, we assume that the propulsion force is a delta function 
-^(x-Xp)Fp. 

We now discuss the dynamics of a single swimmer immersed in a fluid. The kinematic constraints resulting from 
rigid connections between the balls and the location of the force can be implemented mathematically by the equations 

(v„ - v^) • r = (Vp - v^) -T = (3) 

= = = ^, (4) 

where v^, v^, and Vp are the linear velocities of the head, the tail, and the force. The second constraint Q expresses 
the assumption that the balls do not rotate with respect to the axis of the swimmer: co^j, u^, and tu^ are the angular 
velocities of the balls and the axis respectively (hence we use the notation lo without a subscript). The constraints 
(H) and ^ can be thought of as implemented by a rigid rod (connecting x^ , x^ , and Xp ) of negligible thickness and 
mass, hence of negligible drag and inertia. 



4 



From this point on we will use the symbol * to indicate H, T, or P. When the meaning of * can be ambiguous, we 
will explicitly mention the values it is allowed to take. 

The motion of a point x on the surface of the ball B^, (* = T) can be described in two equivalent forms: 

v(x) = v^, + (x — Xp)w or v(x) = + (x — x,)ti;. (5) 

Under different circumstances it is convenient to use one or another of these forms. The connection between two 
forms is given by 

= — Lt X oj, = + Lr X to. (6) 

Note that any values of and uj define a rigid motion of a swimmer. However, v^, v^, and u must satisfy the 
additional rigidity constrains I 
and Bt balls being preserved: 

T-K-vJ = 0. (7) 

The second rigidity constraint is the consistency of the rotation defined by v^^^, v^,, and lj: 

K - V J = -2iT X (8) 

The second rigidity constraint shows that the linear velocities and of the balls B^ and B-^ define the angular 
velocity cj of a dumbbell (up to a rotation around the dumbbell axis t) . 
From Newton's second law of motion 

mv^ =Fh + Ft + Fp, (9) 
/w = Th + Tt. (10) 

Here = x^ is the velocity of the center of mass x^ of the dumbbell, / is the moment of inertia for the dumbbell 
with respect to x^, and Fp = /pT is the reaction force onto the dumbbell from the point force pushing on the fluid 
(modeling the action of flagellum). The forces F^ and F^ are due to the viscous drag exerted onto the head and 
the tail by the fluid; likewise, and are the torques due to the viscous force (the applied force Fp acts along 
the dumbbell's axis and results in zero torque). All forces are applied at the center of mass x^ and all torques are 
calculated with respect to x^. The forces F* and the torques (* = H, T) are given by 

F, = /" cr(u) • n(x)dS'(x), T, = /" (x - x^) x ct(u) • n(x)dS'(x), (11) 

JOB, JdB, 

where the stress tensor (t(u) is defined in terms of the strain rate (symmetrized gradient) -D(u), 

ct(u) = 2^i:>(u) -pi, 2D(u) = Vu+ (Vu)^, 
and n(x) is the unit inward normal to ^B.^, . 



Tpl to define a rigid motion. The first rigidity constraint is the distance between Bh 



B. Discussion of the model 



Many bacteria swim by rotating the flagellum, driven by the torque generating motors located within the bacterial 
membrane. The rotation of the flagellum causes the body of the bacterium to rotate in the opposite direction around 
its axis of symmetry. This rotation exerts a hydrodynamic torque on the fluid. Experiments for which this torque is 
important are those where bacteria are close to one another or to the wall of container. For example it is known [27] 
that a bacterium swimming next to a solid wall will swim in circles in the plane parallel to the wall. The direction of 
the swimming is determined by the chirality of the flagellum. 

On the other hand, for well-separated swimmers the effects of the torque around the axis of symmetry can be 
neglected because the disturbance due to it decays as r~^, faster than the decay of the disturbance due to self- 
propulsion. Ignoring the motor torque simplifies the model, making it readily amenable to analytical treatment while 
still capturing the key features of the experimental observations in the dilute limit. 

Note that our model is in fact a model of a self-propelled swimmer as opposed to a body propelled by an external 
force, such as gravity or the magnetic field. The propulsion of our swimmer is due to the point force that models the 
effective action of flagellum on the fluid. This force is balanced - see equations ( 15 )-( 16 1 ) - by an equal and opposite 
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force of the fluid onto the flagellum transmitted to the body of the swimmer. Therefore the propulsion force is not 
external since all external forces (e.g., gravity) are unbalanced. 

Another issue of concern at very small scales is Brownian motion. How reasonable is it to ignore such motion? The 
rotational diffusion coefficient Drot for an ellipsoid of length I and diameter d (see [551 HH] ) is 

^ 12 fc„r 

J^rot — 



TT Pln{l/d)ri ' 



where k„ is the Boltzmann constant 



fc„ = 1.3806503 X 10" 



kg 



T is temperature, and r/ is the viscosity of the fluid. 

Computing the value of the diffusion coefficient i^rot for Bacillus subtilis (/ w 4 — 5 /im and d w 0.7 — 1.0 /im) in 
water at near-room temperature (T « 300 K and 77 — 0.8 x lO^^N • s • m"^), we obtain 

_ 12 1.3806503 x 10-23.300 m^ • kg • K • m^ _ 2 ,-1 
™* V (5- 10-6)3 .ln(5) •8- 10-3 s2 • K • m3 . N • s ^ ' 

An isolated swimmer in the absence of Brownian motion will swim in a straight line. In the presence of Brownian 
motion, the expected time T{6) to deviate by an angle 6 from a given orientation can be computed by using the 
appropriate flrst passage time as 

T{9) = 



Thus, for "interior" swimmers in the "mirror image" conflguration (see section III C 1 1 the expected time to leave 
the basin of attraction of the "swim in" conflguration due to the described thermal effects becomes comparable to 
the interaction time (the "swim in" time) at distances of 100 /im and larger. At smaller separations we can, therefore, 
ignore thermal effects, at least at the qualitative level. 



C. PDE and ODE models for a swimmer in a fluid 



In the preceeding section we described the dynamics of a swimmer. In this section we derive the PDE and ODE 
models for a collection of swimmers interacting with a ffuid. We consider several neutrally buoyant swimmers (indexed 
by a superscript i — I, . . . , N) immersed in a Newtonian fluid (water) that occupies the domain fl. We are concerned 
with instantaneous velocities of swimmers and fluid due to propulsion forces. The head and the tail balls of the zth 
bacterium are denoted by and B^, respectively. The corresponding coordinates, velocities, forces, and torques 
are labeled accordingly. The swimmers occupy domain CIb — [Ji ^ Bl while the fluid occupies domain flp = \ fJ^. 

Based on the typical swimming velocities and sizes of swimming bacteria, the Reynolds number of the fluid flow 
induced by the motion of the swimmers is usually less than 10"^ (see, e.g., [H]). Then inertia forces on the fluid 
elements are entirely dominated by viscous forces. Ignoring the inertial effects of the fluid in our model is reflected 
by reducing Navier-Stokes equation to Stokes equation, described below. 

In the Stokesian framework, the Stokes drag law is applicable, stating that the viscous drag on each of the balls is 
proportional to the radius R, while the mass of the ball is proportional to R^. Indeed, a neutrally-buoyant swimmer 
has the density of the surrounding fluid p, so its mass is ^R^p. For sufficiently small R and flnite density p and 
viscosity p, the inertial terms mv^ and Ilu in 
of forces and torques: 

= F^+F^+F^, = T^+T^. 

As discussed above, at any time the fluid obeys the steady Stokes equation on the domain il.p determined by the 
instantaneous conflguration of the dumbbells. The fluid is at rest at the outer boundary dfl (the container) and is 
coupled to the dumbbells only through the no-slip boundary conditions on the surface of the swimmers. Therefore, 
given the fixed magnitude fp of the force Fp = rfp and the instantaneous positions and of the balls B^ 
and i?^, their instantaneous velocities = x^ and = x^ are related to the fluid velocity u(x) through the 
incompressible Stokes equation 

|div(ur=''o^'^^'^"""'^^'^'^ ^-^^ = ^\^B, ns= U Bl (12) 

i; *=H,T 



&pO I can therefore be neglected, so that (9|10l reduce to the balance 
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subject to the boundary and balance conditions 



u(x) = 0, X e dn, 
u = X (x-x*)) , on dBl, 

F' + — — F' 
T' + T* = 0, 



container at rest, 
no-slip, * = H,T, 
balance of forces, 
balance of torques. 



(13) 
(14) 
(15) 
(16) 



where fi is the viscosity of the fluid. 

The system (12)-(16l imphcitly defines an ODE initial-value problem for the swimmers. Indeed, given the instan- 

the corresponding velocities and and the fluid velocity fleld u can be 



and orientations 



taneous positions x^ 
simultaneously determined from 



The well-posedness of (12)- (16 



12l-(16) 



solved for u as a linear function of the velocities v 



is established in Appendix |B[ Heuristically, the first equations (|12[)-(|14| can be 

UJ 



^, ~ ■ Doing so eliminates u from the remaining equations where 
and torques T* 



Since (111 is linear in u, F* and also depend 



and w*. Hence, equations ( 15 )-( 16 1 provide a nonhomogeneous linear system for the unknown velocities 



it enters through the definitions (11 1 of forces F 
linearly on 

and tj' in terms of positions xi and given intensity fp of the propulsion forces F^ — t"^ fp. 
Solving this system for the velocities v* and uj^ and using ([6]), we arrive at an ODE system for x' 



i:-v(/p,x^„,x^^,xg, j = \,. 



.N. 



(17) 



In the remainder of the paper we investigate the dynamics of swimmers, deriving appropriate approximate ODEs 
(the fluid is acting only as a mediator of hydrodynamic interactions between the swimmers) . 



III. ASYMPTOTIC REDUCTION OF THE PDE MODEL 



In the dilute limit the problem of determining the drag forces and velocities on the individual balls can be effectively 
approximated by using three classical relations: the Stokes drag law and the basic solutions for the flow due to a point 
force and for the flow due to a moving sphere. We assume that the bacteria are sufficiently long {2L 3> R) and far 
apart (|x* — x^| ^ 2L, i j) so that all the balls and the point of application of propulsion forces are well separated 

(|C-i|,|C + i|~i). 

At a point x the flow due to an isolated propulsion force from the rotation of flagellum — Fp is given by u(x) = 
—G ■ Fp. Here G is the Oseen tensor; see ( |Al I in Appendix A. The velocity of the fluid due to a ball moving with 
translational velocity v in an unbounded fluid domain is given by u(x) = H v. Also, the drag force onto a ball 
moving with velocity v is given by Stokes formula 



70 V, 



7o — GnRfj, 



(see Appendix |A] for the deflnition of G and H) . 
Furthermore, at distances large compared to the radius of the ball R we have 

i?(x) w 7oG(x) for |x| > R. 



(18) 



(19) 



Using (19 1, both the flows due to a point force — <5(x!,)F* and due to translating spheres Bl, — H,T) we can write 



in the same form (20 1 in terms of the Oseen tensor G and the forces F^ (* = H, T, P) exerted by the fluid 

u:(x) = -G(x - x:)Ft, *^H,T, or P. 



(20) 



We use (18 1 to relate the drag force F* (here * = H.T) on the ball B\ (— F^ is the force onto the fluid) to the 

(21) 



velocity of the ball Bl relative to the flow u* , given by (|23 1 

-f:-7o(v: 

Hence, 

„■ 1 



7o 



* = H,T. 



*^ H,T. 



(22) 
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For a given configuration of balls and point forces, remove one ball Bl and replace it by fluid. Then the velocity of 
the the center of the "fluid ball" Bl is denoted by u* . For instance, for the ball Bj^ (i.e., * = H) is given as 



N 



= V u^' (x^ ) + u^' (x' ) + u^' (xM + (x' ) + u' (x* ) 



(23) 



= - E on + G(x;, - x^)F^ + G(x; - x^F^ 

-G(x* - x^ )F^ - G(x'' - X* )F* . 

\ H T ' T ^ H P ' P 

(Similarly we express .) 

Note that is defined only in terms of the forces F^ and positions x*. Hence, if we know all the forces F^, 
equation ( 22 1 together with ( 23 ) will give us all the velocities . In the remaining part of this section we explain 



how to obtain a closed system of 6-/V equations for the unknown forces F^ (* = -ff, T). Since the number of unknown 
components of forces F^ and F^ is 2 x 3 x TV = 67V, the system is, indeed, closed. 

The relation ( 15 1 consists of 3iV equations. Since i? <C i, we have (x — x^) « (x — x^) = ±Lr. Hence the torques 
T* on the balls B\ can be approximated by the moments of the forces on them, so that the torque balance relation 
(16 1 becomes 



r' X (F^^ -F^) = 0. 



(24) 



Since linear operator A^-F := t x F has a ID kernel, the relation (24 1 gives us another 2N equations. To obtain the 
remaining N equations, substitute the equations (22 1 into the rigidity of bacteria equations ([t]), and use expression 
( 23 1 for and . The obtained 6 equations schematically are denoted by 



A (x:,f;,f;) (xI,f;), 



(25) 



where £1 and £2 depend on all positions x^ and forces F* (* H, T, P, i = 1, N). Here Ci and £2 are nonlinear 
in positions x* , due to the terms of the form G(x* — xi), and linear in forces F* . 



We solve ( 25 1 for F^ and F^ in terms of the known F^ and x* . Then and are defined by ( 22 1 and ( 23 1 

v: = V(F^,x^,x;,x^), 

and we obtain an ODE system (x^ = v*): 

x: = v(F^^,x;,x;,x;) 

In particular, for a single swimmer (A^ = 1) simple computations show 



(26) 
(27) 



3R 



IC-il IC+il 



(28) 



where t is the unit vector along the axis of the bacterium. 

Also, when all the forces have been found, one can find the velocity field of the fluid, using the approximation ( 19 1 

«(x) + u;(x)+u;(x)) = 

'g(x - x^J . F»„ + G(x - x^) . F; + G(x - x;) • F' 
This flow for a single bacterium is illustrated in Fig. |2] 



u(x) 



N 

E 



(29) 



N 

E 

4 = 1 



A. Bacterium as two force dipoles 

The above system of 6iV equations can be reduced to a smaller system for N unknowns if one makes the following 
observation. From the balance equations (15), (24 1 and the form of the propulsion force F^ = fpT^ it follows that F^, 
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FIG. 2: (a) Velocity field of the fiuid, computed fi'om (291, around a swimmer with ( — —3. Heuristics: The bacterium is 
moving left to right; hence the head ball pushes the fluid to the right. The force of the flagellum is pushing the fluid to the 
left. Because of incompressibility, the fluid is forced toward the bacterium from top and bottom, (b) Velocity field of the fluid, 
computed from (29 1, around a swimmer with ( — 0. 



F^, and are all coUinear. Indeed, F^ is parallel to r* by definition, and the balance equations (15 1,(24) imply 
that both the sum and the difference of F^ and F^ are collinear with t*. Hence, F^ and F^ themselves are coUinear 
with r\ 

Therefore, for each bacterium there exists a scalar parameter a' such that 



(30) 



Effectively, the parameter a' groups the forces F' in two "force dipole" pairs as illustrated in Fig. |3] Since F^ are 
given, F^ and F^ are completely determined by the scalar a\ 



-Fp=F,+Fh 
< • 



-Ft 
• > 



-Fh 



+ 



<- 



-Ft 



-Fp 



-Ft 



-F„ 



-(l-a)Fp 

+ 

-aFp 



aFp 



(l-a)Fp 



FIG. 3: Schematic representation of a swimmer as two "force dipoles" in terms of and (left) and in terms of Fp and a 
(right). Arrows on the top line are the sum of arrows on the second and third line. 



We obtain a linear system for a* by substituting ( 23 1 into ( 22 1 and using the rigidity constraint ([t]) for v* 

2 



G(x^-x;) + G(x;-x;)- 



7o 



G(x; - x^^) + G(x; - xg + G(x; - x^j - G(x^ - x^j 



G(x; - x^„) + G(x^ - x^^) + G(x; - x^j - G(x; - x^J 



G(x; - x^) + G(x; - x^) - G(x; - x;) - - 



(31) 



Equation (31 1 all the coordinates x*,xi and the directors t%t^ are known; (r*)-^ stands for transpose of the vector 
T*. The system is linear for a*, and the coefficients depend on the positions x^,xi in a nonlinear way. 
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We now return to the ODE, equations (27 1. As noted before, to find the expHcit form for V, one has to find 
and in terms of F^ and x^. Thus, the exphcit form of V is determined by (30l-(31); see Appendix [c| 



B. Perturbation analysis of the ODE for two bacteria 



Notations 

We next consider (27 1 for the particular case N — 2. Specifically, we focus on the dynamics of a coplanar pair of 



bacteria, when t*, i — 1,2, lie in the xy-plane. To describe the relative positions of the two bacteria, we use the 



2nd bacterium 




1st bacterium 



FIG. 4: Coplanar pair of bacteria at a distance a = e ^ apart. Here 4> :— /{(x^ — x;i), i-axis} is the angle between (x^ — xj) 
and the x-axis; 9^ := Z{r^, a;-axis}, 6^ := Z{r^, a;-axis}. 

angles 6^, 0^, and (f> shown in Fig.[4j The angular velocity in this case is characterized by a single scalar (no rotation 
around the axis of the ith bacterium) and is completely determined by the velocities of the balls of the dumbbell: 



uj = u = 



2L 



(32) 



where t]_ is obtained from r' by a 90° in-plane rotation. The velocity of the center of mass is 



2 2 ' 

and, as mentioned before, the dynamics of the swimmer pair is completely determined by (v* , w*). 
Asymptotic expressions for velocities 



(33) 



As explained above, the ODE system (27 1 can be written explicitly (V can be expressed in terms of a%x* , and fp). 
However, this system is too cumbersome for direct analysis. 

But, since we confine ourselves to the dilute limit regime, we can satisfy ourselves with asymptotic expressions 
in terms of the natural small parameter e = a~^, where a is the distance between the centers of the two bacteria: 
a = 1x2 -xj, I > 1. 



terms at the same orders of e = |x? 



We consider the asymptotic expansion (C9) for a". After substituting it into equation (31 1 for a' and equating the 

(see Appendix C2), we obtain 



a 



ao + 0(e^), 



1 / R f R 

where aj, = - f 1 + z(C) - + ^(C) f 



(34) 



where z{C,) is defined by (C19) 



To this end, we consider the asymptotic expansions 



i ^3 



Vie + Vne + v,£ 



(35) 
(36) 
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and substitute them into the LHS of (32l-(33l. In the RHS of (32 1- (33) we express v* using (22)-(23l in terms of the 
drag forces (or, equivalently, in terms of a ) and expand the obtained formulas in e (see Appendix |c| . 
Equating the terms at every order of e results in the following expressions for v*, cj': 

(37) 
(38) 
(39) 
(40) 
(41) 

where j ^ i. 

At the leading order (e*^) ith bacterium swims straight along its axis t* with a constant velocity vq ^ as if there 
were no other bacterium. Also, each bacterium can be viewed as two "force dipoles" ; see Fig. [3j The disturbance 
due to a point force, given by G(x)r, decays as |x|"^. Hence, the disturbance due to a "force dipole," given by 
(G(x + tL) — G(x — tL)) t, which is like a derivative of G(x), decays as |x|^^. Therefore, the first nonzero correction 
in (351 is V2e^. The first nonzero correction to rotational velocity uj'^ appears only at order e"^. Heuristically, this can 
I, since the RHS of (32 1 is like a finite-difference derivative of the vector field, decaying as = 
corrections starting from are in a separable form 



0{e°) : 


^0 


= 0, vj, = vot'; 


0(£i) : 


L0{ 


= 0, vl = 0; 


0{e') : 


^2 


= 0, ^l=A\fp,L,R,fiX,ao)B'{e\9^ 


0{e') : 




= ^^(/p,L,i?,M,C,«o)G'(0\e',0), 


0{e') : 







be seen from ( 32 



Notice that al 



Matifp,L,R,^i,C,ao)Trig{0\9^ 



(42) 



Her e the function Mat(/p, L, R, /i, ao) is determined by the properties (/p, L, R, fi, C, ao) and ao — ct(L, R, C), given 
by (C18l-(C19l, of bacteria and viscosity /i of the fiuid (material properties); the function Trig(6'^, 6*^, 0) depends only 
on the mutual orientations {9^, 6^, (p) of bacteria. 

The separable form (42 1 allows us to study separately two questions: (a) For given material properties, how does 
the dynamics depend on the initial orientations of bacteria? Here we show that swim in or swim off is determined 
by the sign of Tiig{d^ , 9^ , (j)) . (h) For given orientations, how does the dynamics depends on bacterial structure 
(primarily, the position ^ of the propulsion force)? Here we show that swim in or swim off is determined by the sign 
of Mat(/p, L, R, fj., C, ao)- In particular. 



A\fp,L,R,iiX,ao) > 
A'{fp,L,R,iJ,,C,ao) < 



for pushers {C, < 0), 
for pullers (^ > 0) 



(43) 



(see Appendix C 3 ) . 
For bacterium 1 we have 



and 



A\fp,L,R,^i,C,ao) = J^(l-C-2ao) >0, 
B\9\9^,q}) - -2(l + 3cos(2(02-0))) 



cos(0) 
sin((/)) 



C\9\9^,^) 
E\9\9\,p) 



3sin( 



h+'2.92 - 3</>) +2cos(6'i 



■ 4) [5 cos( 

+ cos(0i- 202 + 0)], 
3/p£^(C^-l) ^ 
2567r/i ' 
35sin(26'^ + W"^ - 50) + 5sin(26|i + 9'^ 
+5 sin(26'^ - fl^ - 0) - 4 sin(6'2 -</>)- 
20sin(36'2 - 30) + 3sin(26|i - W"^ + 0). 



30) 



(44) 
(45) 

(46) 

(47) 
(48) 



To obtain the corresponding expressions for bacterium 2, we simply switch the indexes 1,2 and replace with (tt-I- 0) 
since is the angle between (x^ — xj) and the x-axis. Note that from (44), the sign of G* and B*, given by ( 46|45 1 
will give the sign of the first-order corrections to v* and w\ 
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C. Dynamics of two bacteria 



The asymptotic formulas (37l-(41 1 describe the dynamics of a well-separated pair of bacteria. The difficulty with 
interpreting these equations is the number of independent parameters {9^, 6^, (j)), which does not allow having a 
single, comprehensive graph for the trajectories of two bacteria. Therefore, we consider two basic, yet representative, 
configurations (see Fig. |5]) where there is only one free parameter and the remaining parameters are fixed. 

The motivation for the choice of these basic configurations is twofold. First, the evolution of a simple symmetric 
states, such as "mirror image," provides insights into the behavior of the pair of bacteria in the course of collisions. 
Second, these configurations allow for at least qualitative comparison with experimental data (on swim in/swim off 
of a pair of bacteria as shown in [5]). 



2nd bacterium 



1st bacterium 



— 






a=e' 








,<t>='V2 



2nd bacterium 




(a) "Mirror image" configuration 



1st bacterium 

(b) "Parallel" configuration 



FIG. 5: Two basic configurations for the relative position of two bacteria. Configuration 1 is called the "mirror image" 
configuration, because bacteria are symmetric relative to x-axis. Configuration 2 is called the "parallel" configuration, because 
bacteria are parallel to one another. 



1. "Mirror image" 

We first consider the case with the two bacteria positioned symmetrically with respect to the a;-axis (see Fig.|5ja)). 
Because of the symmetry, the positions of the bacteria will remain symmetric relative to the a;-axis at all times. Then 
the factors ( 45|46|48" l in the asymptotic expressions (39)-(41 1 become 

(49) 



and 



B 



B^ 



-9\7r/2) = 



-2(1 ~3cos(26li)) 



,-9\7r/2) = B\9\^9\7r/2) 



E\9\~9\Tr/2) = cos(6'i) [2-56cos(26'i) +6cos(46i^)] 
E^{9\-9\ttI2) = ^E\9\~9\tt/2). 



(50) 



(51) 
(52) 



Analysis of steady states 

If a steady (invariant) configuration of two bacteria exists (determined by 9^), it has to be rotationally steady, 

oj\ei)^0^cu\9i), (53) 

and translationally steady, 

vi(0i) = v2(0i). (54) 
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rotationally steady 
translationally steady 



71 +e9i 




FIG. 6: Figures (a) and (b) correspond to the bacteria for which the position of the propulsion force is given by = —2. 
(a) Schematic illustration of the dynamics of the angle 8^ (angle between the axis of the first bacterium and the x-axis) 
in the "mirror image" configuration. Arrows on the unit circle indicate the direction of change of the angle 8^. The states 
8^ — ±7r/2 (empty circles) are the unstable steady states for pushers ((" < 0) and stable for pullers ((" > 0), from (431. The 



states 6^ = —681 and 8^ — Tv + eOl (solid circles) are the stable steady states for pusher ((^ < 0) and unstable for pullers {(^ > 0) 



from (43 1, (b) The angle 8^ is plotted against the distance |xj, — x^| between bacteria. For 8^ > (bacteria oriented inward) 



the dynamics indicated by arrows shows that bacteria move toward each other (jx^ — x^| decreases) and rotate outward {8^ 
decreases). This action corresponds to the fir st p art (T o < t < Ti) of trajectories in Fig. (TJa). The bold red curve indicates 
the rotationally steady states, obtained from (531 and (40h(41l. The dashed blue curve indicates the translationally steady 
states, obtained from (541 and (37l-(39l. These curves never intersect - no state is rotationally and translationally steady at 
the same time. 



We show below that no value of 6^ satisfies both (53 1 and (54 1 (while each of these conditions can be satisfied 
separately) . 

The trivial rotationally steady states are 0^ — ±7r/2 (respectively, bacteria moving toward or away from one other 
on a vertical line). These configurations are rotationally steady, since the configuration and the PDE ( 12 1 are invariant 
under reflection across the j/z-plane: 



u^{-x,y,z) 



uy{x,y,z) 
u''{x,y,z) 



p{-x,y,z) ^p{x,y,z). 



(55) 



Hence the trajectories and orientations of bacteria will also be invariant under this reflection. Thus, bacteria starting 
with their centers xj, on a vertical line {x = z = 0) and oriented vertically (0^ — ±7r/2) will move vertically on that 
line. Hence, these configurations are, indeed, rotationally steady. 

But, since the distance between bacteria is not preserved, these configurations are not translationally steady. Thus, 
these configurations are not steady. 

Two other rotationally steady angles (both are stable under variations of 6^ for pushers and unstable for pullers, 
due to (43l) are 



6^ ^0-ele + O{e^), and 0^ 



0le + O{e^), 



(56) 



where 



^1 ^ ^ Ayp,L,fi,C^') ^ -3L(C+1)(1-C) ^ Q 



D^ifp,L,fi,0 



2(l-C-2a2) 



The angles ( 56 ) are found by setting 



For these angles 0^, the vertical component of the translational velocity has the form 



Vn sm 



^) + 0(e2) = -|v^| 



^^| 2AU,Yy) s + Q(s^)<0 
Di(/p,L,/i,C) 
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for £ <C 1. Thus bacteria are moving apart and the states are not translationally steady. 

Therefore, there is no steady "mirror image" configuration of bacteria under the assumptions of the model. 



Dependence of the dynamics of bacteria on the position ^ of the propulsion force 
Next, we plot the trajectories of two bacteria in the "mirror image" configuration. 

We choose an initial orientation of bacteria parallel to the z-axis {d^ — 0). The trajectories of the centers of bacteria 
are shown in Fig. [Tj 

We observe that when the propulsion force is applied between the dumbbell balls (—1 < C < 1), bacteria initially 
move apart and rotate inwards (see Fig.[7jc). After time to (when bacteria have rotated sufficiently inwards) they start 
approaching each other and swim in. Eventually, the distance between the bacteria decreases, and the assumptions 
about well-separated bacteria become invalid, so more accurate representations of the drag forces and velocities are 
needed to address evolution of the pair in this state. Remarkably, this behavior is consistent with the experimentally 
observed attraction between two nearby bacteria; see Fig. 1 in [21 . While experiments suggest the existence of long- 
living states of a close pair of bacteria swimming on parallel tracks, it is likely that this state cannot be properly 
captured in the asymptotic far-field approximation for the velocity fields of moving spheres used in our paper. 





(a) (b) (c) 

FIG. 7: Trajectories of initially parallel bacteria in the "mirror image" configuration (radius of dumbbell balls _R = 1). (a) 
External pushers ((" < — 1) at first (To < t < Ti) attract and rotate outward. When rotated sufficiently outward (t > Ti), the 
bacteria swim off. (b) External pullers > 1) swim ofi'. (c) Internal swimmers (|("| < 1) at first (To < t < Ti) repel and rotate 
inward. When rotated sufficiently inward (t > Ti), the bacteria swim in. 



The detailed explanation of this behavior is as follows. Initially, at the leading order (e") the translational motion 
is along the a;-axis (vq = vqt \\ ox). The next correction e^v| is directed outwards; hence, initially the bacteria move 
apart. This can be seen by substituting (37 1- (39 1 and ( 44|50 ) into the expansion (35 1. As the bacteria move apart 
they are rotating inwards due to the s^ul term in (36). The rotation changes the orientation of r* and hence of the 
leading-order translational motion ~ vqt. At t ^ Iq, which solves 

(v^(to)+£^V^(io)) -62 = 0, 



(57) 



the bacteria rotated sufficiently inwards that the terms Vq and £^v| balance each other. After this moment {t > to) 
the contribution of Vq to the motion along oy-axis dominates e^V2. Hence, bacteria start approaching each other 
(swim in). 

We also observe that when the propulsion force is not between the dumbbell balls (|C| > 1), the bacteria swim off 
(see Fig. |7|a) and|7|^b)). These observations emphasize the fact that the dynamics of the pair of bacteria depends 
sensitively on the position of the propulsion force and, consequently, on the shape of the microorganisms and the 
structure of its propulsion. 

Stability of the "mirror image" configuration: 

The "mirror image" configuration is a reduction that allows us to describe the state of the swimmer pair with only 
two parameters (a and 9^). How generic is this subset within the space of all configurations? Appendix [p] address 
this question in some detail and shows that the "outward" configuration {9^ < 0) is stable whereas the "inward" 
configuration (9^ > 0) is unstable, so that nearby configurations in the "general position" tend to approach the 
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"outward mirror image" state but not the "inward mirror image," at least when the interbacterial distance is large: 
a ^ 1. Since the "inward mirror image" is central to our description of asymptotic scattering of swimmers, we briefly 
comment on its validity. We regard this configuration as representative of the general asymptotic dynamics in that 
if a "swim off" (see below) occurs for the interaction of swimmers in the "inward mirror image" position, it will 
certainly occur in the "general position" case. At the same time, if a "swim in" occurs in the "inward mirror image" 
situation (as is shown below for specific choices of the force location), it is likely to occur for the nearby "general 
position" configurations, since the crucial "inward" character of the configuration is robust to perturbations. We plan 
to investigate this matter more closely in the future by considering a wider subspace of swimmer pair configurations. 




FIG. 8: Pictures (a)-(c) show the trajectories for inner swimmers ((" = 0) starting from a perturbed parallel "mirror image" 
configuration. The measure of perturbation 5 is defined by | |D2[ ). The picture (a) corresponds to 5(0) — 0, that is, the 
unperturbed "mirror image." The picture (b) corresponds to 5(0) — 0.01, and the picture (c) corresponds to 5(0) — 0.1. The 
unit of length here is the radius i? = 1 of a baU in the swimmer dumbbeU. 




FIG. 9: Distance between two swimmers, starting from a perturbed parallel "mirror image" configuration distance 200 apart, 
at time T (T ~ 100 seconds) as a function of the initial perturbation 5(0). The unit of length here is the radius i? = 1 of a ball 
in the swimmer dumbbell. 



The "mirror image" configuration is stable under small perturbations in orientations of the swimmers when they 
are rotated outward from one another. The "mirror image" configuration is unstable under small perturbations in 
orientations of the swimmers when they are rotated inward to one another (see Appendix [P]) . 

Nevertheless, the swim in or swim off of swimmers in the "mirror image" configuration (as can be seen from Fig. [s] 
and Fig. ^ is representative of their mutual dynamics, resulting from hydrodynamic interactions. 



2. "Parallel" configuration 

Here we consider a pair of bacteria in a "parallel" configuration (see Fig. [5| : one located ahead of the other parallel 
to one another {6^ = 6'^) and initially parallel to the a;-axis. The parameter measures the angle between (x^ — x^) 
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and the x-axis. As bacteria may not be ahgned with the x-axis for t > 0, we introduce another parameter 
angle between (x^ — and the axis of the first bacterium T^{t). Thus, (j){t) = (j){t) — d^{t). 



the 



The factors (39l-(41 1 in the asymptotic expressions (|35|-(|36|, when written in terms of (p instead of 0, become 



and 



= C\0,0,(j}) =3sm{2(j)) l-5cos(20) 
C^{9\e^,4)) = Ci(0,0,^ + 7r) = 3sin(20) [l - 5cos(20) ^C\0,0,4>) 



B\9\e^,4>) = Bi(O,O,0) = -2(l + 3cos(2^)) 
B2(6»^6'^0) = Bi(O,O,0 + 7r) = 2(l-3cos(2^)) 



cos(0) 
sin(0) 

cos(0) 
sin(0) 



(58) 
(59) 

(60) 
(61) 



--sin(0) 



9 + 2Ocos(2<^) + 35cos(40) 



Since for a general angle cj) term £'^(0, 0, (p) ^ 0, equation (|63]) implies that 



(62) 
(63) 

(64) 



This means that a pair of bacteria in the "parallel" configuration may not remain in the "parallel" configuration at 
some later time. In other words, the "parallel" configuration may not be preserved in time (unlike the "mirror image" 
configuration, which is preserved in time). 

The only angles (j) for which bacteria remain in the "parallel" configuration are (f> — and (j) = tt. These angles 
correspond to a pair of bacteria one following another on the same straight line; we call this the head-to-tail config- 
uration. The difference between (j) = and = tt is only in assigning numbers to bacteria ((/) = tt means that the 
leading bacteria is called the first, while (p — means that the traiUng bacteria is called the first). Next, w.l.o.g. we 
consider the case = 0. 

"Head-to-tail" configuration 

From the top-bottom symmetry, it follows that for (p ~ th e ro t atio nal cor rect i ons at all ord ers vanish. For instance, 
to the order this can be checked by plugging (p = into (40l-(41l using (44 1, (46) and (48 1 

The stability (under variations in (j)) of the "head-to-tail" configuration of bacteria is dete rmined by the sign of the 
leading-order correction terms and in the rotational velocities oj^ and ui"^; see ( |58|59[ ). 

Take 6*^(0) — 9^{0) (bacteria initially aligned with a:-axis) and (p{0) — (j){0) small positive (the second bacterium 
is ahead and slightly above the first one). Then = < 0, which means uj^,uj'^ < 0: the bacteria are rotating 
clockwise. The angle (p = (p — 9^ increases. Similarly, take (j>{0) = 0(0) small negative. Then = > 0, which 
means uj^,uj'^ > 0: the bacteria are rotating counterclockwise and (p — (p — 9^ decreases. 

Therefore, from (43 1, for pushers (^ < 0) the "head-to-tail" configuration is unstable and for pullers (C > 0) it is 
stable. This result is in fact consistent with the simulations of j24j indicating formation of close "head-to-tail" pairs 
of puller dumbbells. 

Dependence of the dynamics of bacteria on the position C, of the propulsion force 

We study the dependence of the dynamics of two swimmers in the "head-to-tail" configuration, depending on the 
position C, of the propulsion force. Since the two swimmers are positioned on the same line (there is no preferred 
direction other than this line), will stay on this line and can either get closer together or get farther apart as they 
move on this line. 

We observe that pushers (^ < 0) always swim off and pullers (^ > 0) always swim in (see Fig. 10). 



D. Quasi-two-dimensional model 



In this section we consider two bacteria swimming in a thin film (quasi-two-dimensional fluid, abbreviated Q2D). 
The interest in studying this case is due to a number of physical experiments (e.g., [I] [21 [51 [^) observing the motion 
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FIG. 10: Dependence of the relative velocities of two swimmers in the "head-to-tail" configuration on position ( of the 
propulsion force. Pushers ((" < 0) swim off; pullers {( > 0) swim in. The asymptotic technique used by us cannot be applied 
to the uncharacterized regions (between dashed lines) of ^ close to ±1. 



of bacteria in a thin film (in particular, in experiments in ^ the thickness of the film was of the same order as the 
thickness of the bacteria) . The thin film allows us to focus a microscope on individual bacteria and track their motion 
with time. 

The modeling in a thin film (of thickness 2h) differs from the above model in the whole space because the boundary 
conditions on the top and bottom of the thin film must be taken into account. While the experiments in [TJ[5] were 
performed with free-standing fiuid film, suggesting free slip boundary conditions on the interfaces, the experiment in 
[5^ indicates formation of thin, solid-like walls on the fiuid-air interfaces due to the byproducts of bacteria metabolism. 
Therefore, in fact, the correct boundary conditions for the in-plane velocities are no-slip. 

Hence, instead of the fundamental solution G(-) of the Stokes equation in the whole space, we use its Q2D analog 
- the Green's function G(-) with no-slip boundary conditions on the horizontal walls {z = ±/i): 



Gix,y,h) = G{x,y,-h) = 0. 
The series expansion for the velocity of the fiuid due to a point force 6{r)ei is obtained in |30j : 



Taking the leading term in this series (as |r| 



u(r) 



u(r) = G(r)ei. 
cxd), we get an approximation 



Uy{r) 



2xy 




(65) 



(66) 



(67) 



where /(z) is a known function (see |53 and Appendix A4| 



Analogously to the 3D approximation (19 1, we want to approximate the fiuid fiow due to a sphere moving (in the 
xy-plane) midway between the walls by 



u(r) « -7oG(r)F, 



The approximation (|68| is valid when i? <C /i 

1 



(68) 

It applies here, because we 



i? < /i < £ < £~ , where e~ 



xi, — x^, I is the distance between 



where F is the drag force on the sphere 
are concerned with the following scaling regime 
the two bacteria. 

The solution procedure is exactly the same as for the 3D fiuid, except that G(r) is replaced by G(r). Using a Q2D 
analog of (29 1, we obtain the velocity field of the Q2D fiuid due to a swimming bacterium (see Fig. 111. 

Note that the asymptotic Green's function G(-) for the Q2D fiuid is qualitatively different from the Green's function 
G(-) in a 3D fiuid. For instance, they have different rates of decay: G(r) ~ \r\~^ and G(r) ~ \r\~'^. In addition, since 
the shear modes in the Q2D geometry decay exponentially with the decay rate determined by the spacing between 
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the wall 2h, (see, e.g., |3T]), only curl-free "pressure modes" decay powerlike survive far away from the origin. But, 
most important, G has negative coupling, efG{e2)ei < 0. This means that by applying force to the Q2D fluid in the 
positive direction along the x-axis some of the fluid will actually be moving in the negative direction (unlike in 3D, 
where the coupling is positive and all fluid moves in the positive direction). In spite of these qualitative differences, 
the velocity of the fluid due to a swimming bacterium in Q2D and 3D fluids have similar structures (compare the 
bold arrows on Figs. |2]and 111. This similarity of the velocity fields suggests that the dynamics of bacteria may also 
be similar for 3D and Q2D fluids. Indeed, we find this to be the case. 




FIG. 11: Velocity field of the Q2D fluid due to a single swimmer: (a) pusher and (b) mid-swimmer. 



Note that the velocity field (67 1 due to the point force is curl- free. Therefore, the velocity fields (a) and (b) are 
also curl-free as superpositions of velocity fields of the form (67 1. It appears that circulation of the velocity fields 
in Figs, [nja) and (b) along closed contours passing through dumbbell balls is not zero, since vector fields point 
counterclockwise along some curves. This situation leads to an apparent contradiction with the Stokes formula. 
However, all such curves pass through a singular point of the vector field in the center of the ball, and the Stokes 
theorem does not apply (compare to classical electrostatics where all field lines pass through point charges). 



Asymptotic expressions for velocities 

Substitute the asymptotic expansion ( 35[)-( 36) into the LHS of (32l-(33). Write the velocities of the balls in the 



powers of e. 



RHS of (32)-(33l in terms of ao; see (CI l-(|C2[) and (C9l. Expand Gj^ in powers of e and solve the equations at like 



0(1) : 




= 0, 






0{s) : 


-I 


= 0, 


-\ 


-0, 


O(e^) : 




-0, 




= 0, 


0{s') : 


L0\ 


= 0, 


-\ 


/pi(l 



47r/i 



-cos(26'2 - 3</)) 
sin(26l2 - 30) 



O(e^) 



3/pL(l - C - 2ao) 



4:77(1 

SUL^ie - 1) 
27r/i 



sin {29^ + 29^ - Aq 
in {29^ + 36*2 



50) 



(69) 
(70) 
(71) 

(72) 
(73) 
(74) 



Next, we analyze the dynamics of two well-separated bacteria in the "mirror image" and "head-to-tail" configura- 
tions (see Fig. |5| in the Q2D fiuid. We observe that the dynamics of bacteria is qualitatively the same as that of a 3D 
fluid. The robustness of the dynamics can be explained by the similarity between the velocity fields due to swimming 
bacteria (compare Figs. |2] and [iT]). 



Dependence of the dynamics of bacteria on the position of the propulsion force 



A. ("Mirror image" configuration, Q 2D fluid) 
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We analyze the dynamics of bacteria depending on the position C of the propulsion force for the "mirror image" 
configuration of bacteria. We observe that (as in 3D, see Fig. ulc)) when the propulsion force is positioned between 



the dumbbell balls (|C| < 1) the bacteria swim in (see Fig. 12 ' 

Also, (as in 3D, see Fig.[7](a) and Fig.|7|b)) when the propulsion force is positioned outside the dumbbell (|C| > 1) 
the bacteria swim off (see Fig. and Fig. ^b)). 
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FIG. 12: Trajectories of two swimmers in a Q2D fluid in the "mirror image" configuration, starting from (To) parallel orientation 
{9i =02=0; radius of dumbbell balls R = 1). Initially in (a) through (c) the bacteria are parallel to each other 6\ = 62 = 0; 
(a) C, < —1, the outer pushers swim off; (b) C > !> the outer pullers swim off; (c) |C| < 1, the inner swimmers swim in. 




FIG. 13: Dependence of the relative velocities of two Q2D swimmers in the "head-to-tail" configuration on which determines 
the position of the propeller. Pushers ((^ < 0) swim off, and pullers (0 < swim in. The asymptotic technique used by us 
cannot be applied to the uncharacterized regions (between dashed lines) of ^ close to ±1. 

B. "Head-to-taiV configuration, Q2D fluid 

For two bacteria in the "head-to-tail" configuration in the Q2D fiuid we observe the same dynamics as for the 3D 



fiuid - pushers (C < 0) swim off and pullers (0 < C) swim in (see Fig. 13 1 



IV. CONCLUSIONS 



In this paper we studied the hydrodynamic interaction between two microscopic swimmers, modeled as self-propelled 
dumbbells, in two distinct settings: three-dimensional and quasi-two-dimensional fluid domains. The interaction in a 
three-dimensional fluid domain models the interaction of swimmers in the bulk (away from the walls of the container) , 
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while the interaction in a quasi-two-dimensional fluid domain models the interaction of swimmers in a thin film. 
Qualitatively, models in both settings produced the same results, thus suggesting that the hydrodynamic interaction 
of a pair of swimmers is robust under the change of geometry of the fluid domain. 

At the same time, the shape of the swimmer, that is, the position (C) of the effective propulsion force, proved to have 
a critical effect on the character of the hydrodynamic interaction of swimmers. In the "mirror-image" configuration 
the dynamics of swimmers differentiates inner (|^| < 1) and outer swimmers {\(\ > 1). Inner swimmers (|^| < 1) in 
the "mirror-image" configuration experience a swim in, approaching each other and perfectly matching the swim in 
experimentally observed in [9 for the rod-shaped bacterium Bacillus subtilis, which has multiple fiagella distributed 
over the cell surface. Unlike inner swimmers, outer swimmers (|C| > 1) in the "mirror-image" configuration experience 
a swim off, due to outward rotation. 

In the "head-to-tail" configuration the dynamics of the swimmers differentiates pushers {( < 0) and pullers (C > 0). 
Pushers (C < 0) in the "head-to-tail" configuration experience a swim off; that is, while they remain oriented along 
the same straight line, the distance between them gradually increases. Unlike pushers, pullers (C > 0) in the "head- 
to-tail" configuration experience a swim in; that is, while they remain on the same straight line, the distance between 
them gradually decreases. Moreover, for pushers the "head-to-tail" configuration is not stable, whereas for pullers 
it is stable. Thus, our model predicts a formation of "head-to-tail" structures by pullers and no such structures for 
pushers. 

The surprising sensitivity of the observed hydrodynamic interaction of swimmers to the fiagellum position (more 
generally to the structure of the propulsion apparatus) and, therefore, to the structure and the shape of the swimmer 
explains the wide range of behaviors exhibited by microorganisms (see, e.g., [32] for a study of a sperm cell with a very 
long fiagellum and f33^ for a study of algae that pull themselves forward with fiagella positioned in the forward part 
of the body) and different models of microscopic swimmers, such as dumbbells |24], squirmers |25j . self-locomoting 
rods |23j . and three-sphere swimmers |34j . 

Further refinements of our model are keenly needed. In particular, our calculations are conducted in the dilute 
limit, where the distance between the swimmers is large compared to their size. However, as we demonstrated, pushers 
have a tendency to converge, thus eventually violating this approximation. Therefore, nontrivial regularizations of 
the interaction at small distances using, possibly, lubrication forces and hard-core repulsion must be included into 
the model in order to obtain agreement with experiments and simulations. Further, at high concentration, deviations 
from the pairwise interaction may also become important, especially because hydrodynamic forces decay very slowly 
in the three-dimensional geometry of the sample. 
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APPENDIX A: BASIC STOKES SOLUTIONS 
1. Point force 

The velocity field due to a point force F in an unbounded fluid domain is 

1 (, xx^ 



u(x)^G(x).F, G(x) = — J. (Al) 

Tensor G (along with a suitable pressure tensor P) solves the Stokes problem with a point force 

div(u) = 

Therefore, it is the fundamental solution to the above problem, given in components by 



' 8^m|x| V ' |x|2 
with the corresponding pressure, a vector, given by (up to an additive constant) 

' ~ 4^|x|3- 

The stress tensor corresponding to G and P is a triadic S: 

3 Xj^x ' X 

^ijk = -PjSik + 2 {Gij,k + Gkjd) = |^|5 ■ 

For more details see [55] , 

2. Swimming ball 

A ball of radius R moving with a constant velocity v through an unbounded fluid domain creates the velocity fleld: 

u(x) = i/(x; i?) V, i7(x;i?)= If [al + /3nn^] , (A2) 

ML 

3^2 : 



where I is the identity matrix and (nn-^) v = (v • n)n is the dyadic product. 
Away from the origin (r 3> R) 

ff(x; R) « 70 G(x), 70 = 67r^i?, 

where 70 is the inverse mobility of the ball, characterizing the applied force necessary to generate a steady translational 
velocity of unit magnitude. 

3. Stokes law for drag 

The drag force from the fluid of viscosity /i on a ball of radius R, moving with a velocity v through unbounded 
fluid is 

F = -70V. (A3) 

More generally, suppose that a ball is added to given an initial background flow u and that under the influence of 
external forces the ball undergoes a steady tranlational motion of the ball with velocity v. The Stokes law for drag 
states that the accompanying drag force F on the ball is proportional to the difference of the velocity of the ball and 
the velocity of the background flow, which would exist at the location of the ball x in its absence: 

F = -7o(v-u(x)). (A4) 
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Since in the Stokes framework the drag on the ball must be balanced by the applied forces on the particle, (A4| 
provides a means of calculating the net applied force that results in a given translation velocity v. 
The Stokes law is an approximation to Faxen's first law [551 : 

F = -7o(v-u(x)) + 7o^V2u(x). (A5) 



If th e ba ckground flow is due to a point force or another translating sphere at x, far from x, then it follows fr om ( Al I 
and (A2| that the gradient is small - ~ jsr^sp- ^""^ '^^^ ^^'^ Stokes law (A4| is a good approximation to (A5). 



4. Q2D Green's function 



Take formula (51) in |30] 



TT^i H \ h) H \ H 



1 J rar/3 



+ Sja^kl3 



o {"^'J^p-he-^y^i^ + o r-^'^ip-h-"^/" 
p p J \p p 



(A6) 



where yi ~ 4.2, and rewrite it in our notations. The point force is applied midway between the walls of the film. 
Replace h — ^H; here H is thickness of the film (replace by h). Assume k — 1, that is force is applied along ei. Here 
p is the radius vector from point force (replace by r = |r|). Replace x^, by z. 
Performing the above changes, we obtain 



Ap^irp V 2h 



1 . raVp 

2 p'^ 



-pyi/{2h) 



O 



p p J \ p p J 



Note that only the first term in ( A7 ) does not decay exponentially in p: 

3 





1 . rarfj 


) p2 





47t/' 2h) p2 
AiTfi^ V 2h) p2 



P ^ 



ryri 

n'2 



Rewriting u'^ in components, we have 



Ul 

U3 



Anil 



\ 2hJ «2 



p2 2' 

r2ri 1 r 

- 2^21 

rsri 1 r 

-Tfi 9O3I 



Note that p2 



y . Hence, (A8) takes the form 



2h) / 



4:77(1 



? 9 

— y^ 
2xy 
2xz 



V 2h) fl2 



xz 



Since \z\ < h <^ 1, one has an approximation for the Green's function in Q2D: 

r.2 



G(r) « f{z) 



X' — ?/2 2xy 
2xy — x^ 




where f(z) = z 

^ ' 8tt(i 



Here f{z) satisfies no-slip boundary conditions: f{z = 0) = and f{z = 2h) = 0. 



(A7) 



(A8) 



(A9) 



(AlO) 
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APPENDIX B: EXISTENCE AND UNIQUENESS 



In this section we prove the existence and uniqueness of solutions to the system (( 12 )-( 16 1) under different assump- 
tions on the regularity of the propulsion forces and the size of the container. The most restrictive case of smooth 
forces and a bounded container yields the clearest proof that is essentially classical, but contains a few novel features. 
The other cases refine the argument in the case of point forces and an unbounded container. 

For the sake of clarity we consider the case of a single swimmer. The extension to a multiswimmer system is 
straight-forward. We always assume that the swimmers do not overlap, which in particular excludes the overlap of 
any propeller with any head or tail. 

1. Regular case 

To elucidate the main issues in the existence and uniqueness proof we initially consider the case of regular problem 
data. Assuming that the propeller force density is smooth and the container f2 is bounded, the argument is an 
adaptation of the classical techniques based on coercivity to deduce the existence of weak solutions, followed by an 
application of elliptic regularity results. The novel feature is the presence of a boundary condititon on the forces and 



torques ((15l-(16)) and the consequent use of Korn's inequality in place of the usual Poincare's inequality. 



a. Space of admissible flows 



In order to obtain an appropriate weak formulation of the problem that incorporates the balance conditions ((15)- 
p6|), we have to circumscribe the space in which the solutions will be sought. Consider the following space 

V = {ue (P(nF))^| u(x) = v" + tj" X (x-x^) for xe uas^, u|ao = 0,v",tj" e M^}, (Bl) 

where T) (p,p)) is, as usual, the calss of restrictions to {^Ip) of C5"(M'^) - smooth functions on with compact 
support. 

An equivalent definition of V is through constrains 

u(x) = X (x-x„) nedB^, u(x) = v^-Hw" X (x-x^) x e (B2) 

where v" , v" e are given in terms of v" and by 

v;^=v^-|-..x(x„-xj, v^=v^-Hc.x(x,-xJ, v^ = ^^^±^. (B3) 

or after substituting (x^ — x^. ) = Lt and (x^ — x^) = —Lt 

v" = v" + Lcj X T, v" = v" - X T, 

H C ' T C ' 

Now define V = V a.s the closure of V in H^{ilp). 

It is important to note that V is not empty: the boundary data of the form 



^0, xedn 

can be continued to W(x) defined on all of fii? = \ f^s- To show this we can apply the standard theorem found. 



for instance, in [55], as soon as we show that the boundary data of the form (B4| satisfy the compatibility condition 



/ w(x) • n(x)(ix = — / w(x) • n(x)dx + \J / w(x) • n(x)(ix, 

Janp J an ^ ^ Jasi 

where for consistency with the rest of the paper we denoted by n is the normal pointing into the fluid (i.e., the negative 
of the outward normal). The first integral above is clearly 0, while the second integral vanishes by the divergence 
theorem. Indeed, on the boundary dBjjU dB^ of a fixed i — th dumbbell the field is w(x) = v* + x (x — x^), 
which clearly can be continued into the interior of the dumbbell where it is divergence-free. 

Note that v"^, uj'^ as well as and are uniquely defined for any w in V. For u € V that solves ((12l-(16l) (a 
strong solution) v" and are the linear and angular velocities v and uj entering into ([t]) . Taking a fixed extension 
W for each distinct boundary field defined by a pair (v'*^, uj'^) G M.^ (BM.^ = M.^ we can inn fact showt that V is 
isomorphic to the direct sum ® Vq. Here Vq is the closure in Hq{^Ip) of divergence-free functions from ^{ilp). 
The isomorphism is established by sending w G to ((v'*^, cu'^), w — W) e ® Vb. 
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b. Weak form 



Multiplying (B5l by w e V 



In this section we derive the weak form of ((12 )-( 16 1) on V. To this end, assume there is a strong solution u G V to 
((12l-(16l). Starting with the usual formulation of the Stokes' equation in terms of the Laplacian to derive the weak 
form eventuUy leads to bilinear boundary terms, which are har d to bound from below to show coercivity. Instead, 
adding = /iV(divu) = /idiv(Vu)^ to the first equation in (12 1 we obtain an equivalent formulation in terms of the 
symmetrized gradient D{u): 

2fi div {D{u)) = Vp - F. (B5) 

(2^div (D(u)) + F - Vp) • w 
and integrating by parts yields 

—2^ / -D(u) : Vw da; + 2/1 / n • D{u) ■wdx + + F ■ w dx + pdiv(w) dx — pw ■ n dx = 0. 
Using the incompressibility condition, the symmetry of D(u) and rearranging terms we obtain 

4/i / ^(u) : D{-w) dx^ F -w dx+ n - [2^iD{u) - pi] ■ w dx. (B6) 



At this point the idea is to extend (B6l to V thereby obtaining a weak form of (([12|)-(|16|). However, the presence 
of a bilinear boundary integral term will resist easy lower bounds needed to show the coercivity of the problem. This 
difficulty is circumvented by using the balance conditions ((15l-(16l) to eliminate u from the integral, converting it 



into a linear functional of w. To this end, note that the boundary integral contains a product of w with the boundary 
tractions n • ct(u) generated flow u. 



Rewrite the boundary integral in (B6 ) using the boundary conditions (Bl ) and (B2 1 on w featuring in the definition 
of V: 



/ n • (t(u) • w dx 

Jdnp 



dB^UdB^ 



n • (t(u) • (v^ + tj^ X (x — Xp)) dx 



[ n-(7(u) • (v^+t^^ X (x„ -X^) + CJ^ X (x-x„)) + 
h / n-a(u) • (v^ + tj^ X (x^ -x^) + cj^ x (x-x^)) dx = 



F ■ v"" 



Rewrite (Bll ) as 



(F,+FJ 



(B7) 

(B8) 

(B9) 

(BIO) 
(Bll) 

(B12) 



where F* and T*, * = H,T are the forces and torques associated with u as defined the previous section. We will 
show now that the first term defines a continuous linear functional on w, and the last two terms vanish due to the 
conditions on F*,T* and vj^. Indeed, since r • (v^ — v^) = 0, (i.e., (v^ — v^) _L r), we have 

(v--v-) = rx(rx(v:-v-)). 



Hence, from the relation (|B3|) and the balance of torques (|16| we have 

w _ 



(T,+TJ.^--(F„ 



Lt X + {T^+T^)-uj^ 



[(F,-FjxLr + T„+T,] -c^-^O. 
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Finally, from the balance of forces (T5| we have 



(F„+F,) = - 



F = -v"^ F 



which defines a continuous linear functional of w in terms of the fixed total force Fp. 
Thus the solution of the Stokes equation ( 12 1 satisfies the following variational problem 



where 



4a(u, w) = fe(w) Vw e V, 



a(u, w) := /i / ^(u) : -D(v) dx 



6(w) 



/ (F • w - F^ • v^) 



(B13) 

(B14) 

(B15) 
(B16) 



We want to emphasize that it was the use of the symmetrized gradient D{u.) in place of the usual gradient that lead 
to a boundary integral in terms of tractions, which, apart from having a clear physical meaning, enabled elimination 
of u with the help of the force balance conditions. 



The minimization problem corresponding to the variational problem (B14l is 



where the energy functional is 



mini^ful, 



E[u] = 2a(u,u) - b{u). 



(B17) 



(B18) 



The quadratic term a(u, u) is the usual viscous dissipation rate and the linear term 6(u) represents the work of the 
forces in the fluid and on its boundary. The interpretation of 6(u) as the work done by the forces becomes clearer 
once we rewrite it as: 



6(u) 



F • V da; 



F„ • 

H H 



c. Existence, uniqueness and regularity 



The existence and uniqueness of minimizers of (B17l is proved in a standard way provided that the coercivity of 



the bilinear form a(-, •) can be shown. The coercivity proof, using Korn's inequality, is essentially contained in [37] as 
we now explain. 

Theorem B.l. The bilinear form a(-, •) is coercive on V with respect to the norm \\ ■ ||i induced from H^{flp). In 
particular, a(-,-) defines an equivalent inner product on V. 



Proof. Coercivity of a(-, •) relies in an essential way on Korn's inequality: 

a(u, u) 



|U||2 >C||U"2 



1: 



(B19) 



for some c > (here || • || denotes the L2 norm). The proof of (B19l found in [37] applies to the case for any subspace 
U C H^^flp)) consisting of functions with a zero trace on a part of the boundary with nonzero two-dimensional 
measure. This applies to V as its elements vanish on dfl - the no-slip boundary conditions on the outer boundary 
of flp- In particular, a(-, •) is nondegenerate, since the nontrivial kernel of D(u), consisting of the rigid motions 
u(x) = Uq + X X, is excluded from V due these boundary conditions. The result (B19l is nontrivial, since the 
left-hand side contains only symmetric combinations of the derivatives of u. 

The coercivity proof is completed by showing the existence of the following bound: 



a(u, u) > d\\u\ 



(B20) 



for some d > 0. This replaces Poincare's inequality in the case of the symmetrized gradient. It can be proved for V 
as is done in [37], using the compactness of the embedding V ^ i2(f^F)- This embedding is induced from the usual 
compact embedding H^{ilp) ^ L2{i^F), since V, being a closed subspace of H^{ilp) is also weakly closed (see, e.g., 

m)- □ 
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With the coercivity of a(-, •) proved, the existence of minimizers for (B17l can be proved by standard techniques. 
Since each minimizer satisfies ( |B14 1, the difference of any two of them is a(-, •)-orthogonal to a dense subset of V, 
hence is zero, which proves uniqueness. 



Finally, the unique field u that solves (B17 1 is a weak solution of the Stokes equation on a regular bounded domain. 

]) there exists a unique pressure field p € L2{^f), which 
Since the boundary oiflp and the righ-hand side of (([l2|)-( 14 1 ) 



Therefore, once again by the standard theory (e.g. 
together with u satisfies the a priori L2 estimates [36 
are smooth, these estamates imply that (u, p) are smooth too. By reversing the steps leading to the weak formulation 
(B14), we now see that (u, p) form a strong solution of ((ll2])-( 16 1) 



Point forces 



The limit of point forces (5 ^ is useful because it simplifies many concrete calculations in the asymptotic analysis 
of the model. Heuristically, smooth forces can be replace by point forces because the fluid velocity in i3p is ill-defined 
anyway, being a simplified representation of a complicated periodic action of the flaggelum. Therefore, the precise 
value of the velocity and pressure near the propeller Bp do not matter and can be left undefined. At the same time, 
away from the propeller Bp, both a smooth force density in 'D(Bp) and a point force density produce comparable 
results, as will be shown below. 

The case of point forces, however, does not fit into the existence proof of the previous subsection because F is no 
longer in H^{ilF). In this case, however, there is still a unique solution to ((12l-(16)), regular away from xp, which 
can be constructed with the help of the Green's function for Up. 



Assume for the moment that there exists a unique Uq, the solution to (12 1 with homogeneous boundary conditions 



u(x) = 0, xe^r^F. 



(B21) 



Then, the existence and uniqueness of solution u to ((|12[)-( 16 1) is equivalent to the existence and uniqueness of 
Ui = u — Uq, the solution to ((12l-(16l) with F = and the balance conditions modified to account the forces and 
torques due to the point force fiow Uq. Specifically, the balance conditions (( 15 1-( 16 1) are replaced with the following: 



FF* 

T 

rjri 



-0, 
-0, 



balance of forces, 
balance of torques, 



(B22) 
(B23) 



where Fg ^ and Fq ^ are the hydrodynamics forces and torques on the balls Bl, * = H,T due to a constant flow 



Uq and computed using Uq in place of u in (11). Now the method of the previous subsection applies to the system 



satisfied by Ui with the only modification: the linear functional defined on V by (BI6I is replaced by -Li(-) 



61 (w) := - (F, + F„,„ + F„,J . V- - {Tl^ + T^J ■ 
To complete the proof, it remains to show the existence of Uq, which is done in |B 3 



(B24) 



Green's function for Qp 



Here we briefiy indicate how to show the existence of the Green's function (the Green's tensor for fip). The result 
is well-known and we include it for the sake of completeness. The Green's velocity tensor Q and the corresponding 
pressure tensor V are analogous to the (G, P)-pair of (|A 1| in that they solve 



(B25) 



r fiAGi- ~ xo) = VP(- - xo) - <5(- - xo)I 
\ div(u)(--xo) =0 



with x and xq in f2p and subject to the homogeneous boundary conditions on fip. Once the existence of (Cf, V) has 
been shown, the existence of Uq used in|B]is trivially established: 

Uo(x)-/p^g(x-x»p).r'' 

i 

The sought for Green's tensors are constructed by canceling the boundary values of G and P on dilp as follows: 



g = G-G, V = P-P, 
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where G and P solve ( B25 ) with the zero right-hand side and the boundary conditions 

G(- - xo)|anF = G{- - xo)|sf2F, P{- - xo)|c 



P{- - xo)|anF- 



The existence and uniqueness of (G, P) is easily established by standard methods (e.g., [36]) both in the case of a 
bounded and the exterior Q^p = ^^\^b- The only requirement in the bounded case is the compatibility condition 



G(x - xo)dS'(x) • n(x) = 0. 



This equality easily follows from the divergence theorem applied to G, whose divergence is zero in Li(yip), as the 
following simple calculation shows (summation on j implied and |xp = XjXj): 



X 1 X J 



1 



3 Iy|5 



|X|' 
|G, 



< G 



1 



1 

G = const. 



APPENDIX C: ASYMPTOTIC FORMULAS 

The velocities of the head and tail balls in terms of {a-'} and xi are 

-pG(x;-x;) + G(x;-x;) + l/ 
L 7o - 



Note that 



and 



Hence, 



Wt = E [(1 - - <) + - ^i) - ^(x; - x^j 

+ [(1 - a^)G(x; - xj^) + G(x^ - x^) + -I 
L 7o ■ 

G(-x) = G(x) Vx e IR3 
G(-x) = G(x) Vx e 
G(x;-x'^) = G(x^-x^) = G(2Lr). 



Using the ( C1|C2 1 and the relation 



V* = -(v' +vM, 

C ^2 H ^ 



2L 



(v* - vM X 



we obtain 



JP 



[(1 - a^) {G(x; - x^J + G(x; - x^} 



W{G(x;-x^j + G(x;~x^j} 
-{G(x;-x^j + G(x;-x^J} 



G(2LrO + G(x'^ - x^) + G(x; - x^) + -/ 

70 ■ 



(CI) 



(C2) 

(C3) 
(C4) 
(C5) 



(C6) 
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2L 



T X 



^ [(1 - a^) {G(x; - x^J - G(x; - x^J} 



W{G(x^„-x^J~G(x;-x^J} 
-{G(x;-xg-G(x;-x^J} 



+T* X 



(2a» - 1)G(2LtO + G(x'„ - x^) - G(x; - x^) 



(C7) 



1. Expansion of G{-) 

The expansion of in orders of e is due to the expansions of the Green's function G, e.g. 

G(x* - xn ^ G(x^ - x' + x' - x^' + x^' - xn = 

^ H P ' ^ H C C C C P' 



= G (x' - xM + (x^ - xn + fx' - xM 

\^ H C ' ^ C. P ' ^ C C ' 



(C8) 



Here the quantities (x^ — x'^) and (x^ — x-^) measure the distances in the same bacteria, hence they do not depend 
on e = |x* — x^l^-'. The only quantity that depends on e is (x* — x^). 



2. Asymptotic expansion for a' in powers of e = ]x 



c c 1 



Consider the system (31 1 for a* 



Consider the asymptotic expansion 



a ~ a + eai + e an 



(C9) 



and substitute it into (31 ), where all terms are expanded in asymptotic series in e = |x^ — x^ |^^. Note that the terms 
like G(x^ — x^) do not depend on e; hence, their expansion will have only e° order term (itself). On the other hand, 
the terms like G(x* — x^ ) are of order and will not have order e° terms. 
Thus, at the order the equation ([3ll becomes 



G(x;-x^) + G(x^-x'^)- 



7o 



G(x; - xj^) + G(x';, - x^) - G(x; - x;) 



Express the arguments of G(-) in terms of t*, L, and C,: 

x' -X* ^2Lt\ x^ -x' ^{1~C)Lt^ 
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(-1 - C)Lt^ 



and substitute back into (|C10| to get 

G{2Lt') + G{~2Lt') 



2 
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G(-2LtO + G((l - QLt') - G((-l - C)Lt') 
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Using the definition of 



and the properties ( |C20||C21| , simplify ( |C12[ ): 



1 1 1 

2L ^ 2L " ^ 
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1 

47r/i 



1 



1 



1 



2L 



1 



\l-C\L \l + C\L 2R 



(CIO) 



(Cll) 



(C12) 



(C13) 



(C14) 



Multiply through by iirfiRL: 



-Oin 



R 
L 



1 R R 



R 



Pull out i from the RHS, and solve for 



2 2L |l-C|i |l + C|i 



i-CI li+CI 



Denote 



and perform the expansion of (C16) in terms of ^ 



1 



an 



1 + + 



where 



<■ 4 



z = z(C) := 2 



li + CI li-CI 



4C 



-4C 



ifC<-l, 
if 1< C < 1, 
if 1< C- 



Lemma C.l (Properties of G(-)). Note the following two properties of G{-): 
1. LetqeR and t E M^. |t| = 1. Then 



G{qT) = i^G(r). 



2. Let T £ R^, \t\ = 1. Then 



t^G{t)t = 



47r/z 



Proof. Property 1 follows simply from the definition of G(-): 



G{qT) 

Property 2 follows by simple substitution: 



87r/i|(j| 



T^G{T)r = -^r^ (I + rr^) r = -^(1 + I) = ^ ■ 



3. Sign of A'{fp,L,R,fi,ao) 
Lemma C.2 (Properties of (1 — ^ — 2ao)). Assume fp>Oin the expression 

SZTTfl 



where ao is given by (C18 C19). For pushers (( < 0) 

A'{fp,L,R,fi,ao) > 0. 

For pullers (C > 0) 

A\fp,L,R,ii,ao) < 0. 
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Proof. Use the formula (C18) for ofg to rewrite 



1 - C - 2a^ = 1 - C - 1 + 2^ + 



= -c+ 



(C27) 



Con sider three case: C<^li~l<C<lil<C- For each of these cases use the formula (C19 1 for z(C) to evaluate 
( [C27| . 

Case C, < —I: Here z(C) = 72^, and 



-c 



1-e 



4 ^ 



(C28) 



This expression is always positive when (— C — 1) ~ 1 as^^O. 
Case — 1 < C < 1: Here z(C) = ^^Izi! and 



1-e 



4C i_ 



= -c ( 1 . .2 



4 ^ 



(C29) 



This expression changes sign from positive to negative only as C passes through 0. So it is positive when C < and 
negative when ^ > 0. 

Case 1 < Here z{() = ^r^^ and 



c + 



1-e 



c- 



4c e 



This expression is always negative when (C — 1)^1 as ^^0. 



(C30) 

□ 



APPENDIX D: STABILITY OF THE "MIRROR IMAGE" CONFIGURATION 



Before analyzing the stability of the "mirror image" configuration we determine the quantity, which does not depend 
of the orientation of the oS-axis (i.e., it does not depend on the choice of observer), that characterizes how close a 
given configuration is to a "mirror image" configuration. Then we perturb this parameter by a small amount and 
check whether this parameter is decreasing. If this parameter is decreasing, the configuration is stable; otherwise it 
is unstable. 

Note that for two swimmer, in the "mirror image" configuration 



2(1)- (0^ + 0^) = 27rn, 



(Dl) 



for any choice of the oa;-axis. Moreover, if equation (Dl) holds, then two swimmers are in the "mirror image" 
configuration. Therefore, the quantity 



■ + 24> - {0^ + 0^) 



(D2) 



can be viewed as a measure of deviation from the "mirror image" configuration. 

To perform the stability analysis, we perturb the "mirror image" configuration. That is, we choose |i5(0)| > small, 
and we check whether \S{t)\ decreases with time. If i5'(0) has the opposite sign to (5(0), then \S{t)\ decreases with time 
locally, and the configuration is stable; otherwise it is not stable. 

We have 



S' = 2(1)' - (6*1' + 0^') = 2(1)' - {lo^ + u;"^). 



(D3) 



The expressions for ut^ and can be found from ( |40|44|46[ ) . The expression for 0' can be found simply as a projection 
of the translational velocity difference (v^ — ) onto the unit circle 



— sin(0) 
cos{(l>) 



evo (-sin(0)(cos(6|2) - cos{0^)) + cos(0)(sin(6|2) - sin(6l^))) + ©(e^). 
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where we used 

= ^;ot' + 0(£2), i = l,2. 

Without loss of generality, choose the a;-axis so that tt + = 0, that is, (f) = —n/2, which means the second 
swimmer is directly below the first swimmer. Then 



6' = svo (cos(02) _ cos(0i)) + 0(e3) = -2euosin sin + 0(6^) 

= -2ewosin (^6^ - sin + 0{e^). 



(D4) 



Thus, for swimmers rotated outwaxd (0 < 0^ < tt) the "mirror image" configuration is stable and for swimmers 
rotated inward (0 > > —it) the "mirror image" configuration is unstable under small perturbations. The results of 
the stability analysis are not affected by the type of the swimmer; they are the same for all values of 
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